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Abstract 

Quantum entanglement plays a central role in quantum information processing. A main objective 
of the theory of quantum entanglement is to classify different types of entanglement according 
to their inter-convertibility through manipulations that do not require quantum communication. 
While bipartite entanglement is well understood in this framework, the classification of entangle- 
ments among three or more subsystems is inherently much more difficult. 

In this paper, we study pure state entanglement in systems of dimension 2 (g> m tg) n. Two 
states are considered equivalent if they can be reversibly converted from one to the other with a 
nonzero probability using only local quantum resources and classical communication (SLOCC). 
We introduce a connection between entanglement manipulations in these systems and the well- 
studied theory of matrix pencils. All previous attempts to study general SLOCC equivalence 
in such systems have relied on somewhat contrived techniques which fail to reveal the elegant 
structure of the problem that can be seen from the matrix pencil approach. Based on this method, 
we report the first polynomial-time algorithm for deciding when two 2 ® m (g> n states are SLOCC 
equivalent. Besides recovering the previously known 26 distinct SLOCC equivalence classes in 
2 ® 3 (g) n systems, we also determine the hierarchy between these classes. 



1 This work was supported in part by the National Science Foundation of the United States under Awards 0347078 
and 0622033. 



1 Introduction 



The feature that most distinguishes multipartite quantum systems from their classical counterpart 
is their ability to be in so-called entangled states. Not only does quantum entanglement enable 
apparent "spooky action at a distance" between separated systems |TT|. but it also has the poten- 
tial to fundamentally change and dramatically improve the current information processing and 
cryptographic technologies |20|. It becomes natural then to treat entanglement in a multipartite 
system as a information processing resource that takes on different forms as the system realizes 
different states. Much effort has been devoted to formally quantifying the amount of entangle- 
ment contained in a given quantum state with the motivating intuition being that states having 
more entanglement possess a greater degree of computational and communication power than 
those with a lesser amount. 

Under this interpretation, one may reasonably argue that a state 1 4>) possesses no less amount 
of entanglement than another state \ip) of the same system if the system can be converted from 
\4>) to "free of charge/' or without needing any further entanglement to facilitate the trans- 
formation. The physical operations implementing such transformations is the celebrated class of 
Local Operations with Classical Communication (LOCC), which, as described by its name, con- 
sists of all operations in which each multipartite subsystem is manipulated locally but perhaps in 
a manner globally correlated through classical communication. Thus LOCC has become a major 
framework for quantifying and classifying entanglement (see, e.g., the recent surveys on quantum 
entanglement by Horodecki et al. [17J and Gtihne and Toth [ 13 1). 

Following the notation of Bennett et al. 0, we write \ip) < L occ \4>) if \4>) can De converted to 
\ip) through LOCC. When can a state be converted to another? What's the largest ratio one can 
convert multiple copies of a state to multiple copies of another? When is there a state maximum in 
the sense that all other states in the systems can be obtained from this state? Those are examples of 
the many natural questions that arise. We often have answers for the bipartite case. For example, 
a remarkable theorem by Nielsen f23fl states that for bipartite states, < L occ |</>) if and only if the 
spectrum of \4>)'s reduced density operator (i.e. the Schmidt numbers) is majorized by that of 
Another important observation, made by Bennett et al. [2], is that (for any multipartite systems) 
if two states are equivalent under LOCC, they are related by a local unitary (LU) transformation. 
Thus LOCC equivalence classes are simply the orbits of local unitary operations. Such a partition- 
ing is too fine for most interests: even in the two-qubit case, there exists an infinite number of LU 
equivalence classes. 

However, if the required success probability of both the forward and reverse transformations 
is reduced to be simply nonzero, a much coarser partitioning is achieved. General LOCC trans- 
formations occurring with a nonzero probability are called stochastic (SLOCC) and denoted by 
|^) <slocc \ 4>) if the transformation is from \<j>) to |^). It turns out that SLOCC equivalence classes 
are precisely the orbits under local invertible linear transformations [10 j . Similar to the situation 
with LOCC, bipartite entanglement is well understood under SLOCC. Indeed, for bipartite pure 
states, |V>) <slocc \4>) if an d only if the rank of the reduced density operator of (i.e. the Schmidt 
rank) is no larger than that of \<fi). Thus two states are SLOCC equivalent if and only if they have 
the same Schmidt rank. The optimal success probability can also be computed easily from the 
Schmidt numbers |25|. 
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In contrast, entanglement among three or more parties behaves fundamentally different from 
bipartite entanglement. For example, while there is a maximum SLOCC equivalence class for 
bipartite systems of any dimension, there exists two maximal equivalence classes for the simplest 
tripartite system of 3 qubits fflu| . In contrast to Nielsen's Theorem and the rank criterion for 
bipartite SLOCC conversion, deciding SLOCC convertibility in general encodes many difficult 
computational problems. For the general tripartite conversions, the problem is NP-hard (observed 
in (7j using a NP-hardness result by Hastad |[T6H on computing tensor rank). For converting a 
tripartite state to a bipartite state, it is equivalent to the important problem of Polynomial Identity 
Testing [8j- For certain tripartite asymptotic conversion the optimal conversion ratio is precisely 
the exponent of matrix multiplication J7|. In view of those results, a simple criterion or an efficient 
algorithm for checking SLOCC convertibility or equivalence could be found only for systems of 
restricted dimensions. 

This article studies the SLOCC equivalence classes of tripartite pure states in systems of di- 
mensions 2 (g) m (g> n. Dur et al. presented the first major result in the study of multipartite SLOCC 
equivalence classes by showing there to be six different classes in 2 ® 2 (g) 2 systems |T0| . Their 
work was extended to four qubit systems by Verstraete and co-workers in which already an in- 
finite number of equivalence classes exist fl24"| . For an arbitrary number of subsystems, Miyake 
has shown how multidimensional determinate theory can be used to obtain general properties 
and results concerning SLOCC equivalence [21 j. Specific to tripartite 2g)20n systems, Miyake 
and Verstraete have also completely characterized the equivalence class hierarchy and found that 
for n > 4 exactly nine different classes exist ||22| . Using the method of successive Schmidt de- 
compositions, Cornelio and Piza obtained partial results concerning the equivalence classes in 
2 <g> m n systems 0. Chen et. al completed the finite orbit picture by enumerating all 26 equiva- 
lence classes in 2 (g> 3 (g> 6 systems, and showed that for the 3 (g> 3 ® 3 and 2 ® 4 <g> 4 systems (and all 
systems of higher dimensions) there are infinite number of SLOCC equivalence classes. They used 
a technique called "the range criterion" (6), which states that two states are SLOCC equivalent if 
and only if the ranks of the reduced density operators are identical and their supports are related 
by local invertible linear operations. While these results are quite interesting, the tools used to 
obtain them appear rather ad hoc and neither the criterion in Ref |6] nor any previous technique 
provides an efficient algorithm (or any algorithm at all) for determining SLOCC equivalence. The 
non-invertible hierarchy among the 26 classes has also remained an open problem. 

The main insight of this article is that the theory of matrix pencils is the perfect tool for ana- 
lyzing SLOCC equivalence in 2 ® m (g> n systems. For two matrices A,B£ C mxn , the linear matrix 
polynomial XA + \xB is called a matrix pencil. Two pencils XA + fiB and XA' + fiB' are equivalent 
if there exists invertible P and Q such that P(XA + /j,B)Q = XA' + fiB'. A fundamental result 
is the existence of a canonical form, discovered by Kronecker (see, e.g., Gantmacher UHJ). The 
theory of matrix pencils remained an important subject of study for its applications in control and 
systems theory. An example is the computation of the generalized eigenvalues (see, e.g., Section 
7.7 of |15|)j^] The efficient computation of the Kronecker canonical forms, other canonical forms, 
and related problems is still an active field of research (see, e.g. |[lj and following articles). 

The connection with our problem is that each state ina2§m§n space can be represented as 

2 For two matrices A and B, a vector x and a constant A, if Ax = \Bx, A is a generalized eigenvalue of (^4, B) and A 
the associated eigenvector. The set of generalized eigenvalues are precisely {A//i : det(Ayl + fiB) = 0}. 
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a matrix pencil (see Section[3]for details). The local operations on the second and third subsystem 
brings the corresponding pencil to an equivalent one. While actions on the first subsystem may 
bring the pencil to an inequivalent one, we show that if two states are SLOCC equivalent, there 
are only a small number of choices for operations on the first subsystem that make the result- 
ing pencils equivalent. As a consequence, we derive the first efficient algorithm (in fact, the first 
algorithm) for determining SLOCC equivalence in general 2 ®m®n systems. For the systems hav- 
ing a finite equivalence classes, we derive the equivalence orbits (which was known before) and 
represent them using Kronecker canonical forms. We also determine all possible non-invertible 
transformations among those equivalence classes. 

The rest of this article begins with a brief introduction to some main results in matrix pencil 
theory. We then develop the relationship between tripartite pure states and matrix pencils which 
allows us to derive necessary and sufficient conditions for the SLOCC convertibility of 2 (g) m n 
states. From there, we develop a SLOCC hierarchical schematic of all tripartite systems possessing 
a finite number of SLOCC orbits. The article closes with some brief concluding remarks. 

2 Matrix Pencils 

The theory of matrix pencils was first developed by Kronecker over a century ago. A completely 
thorough treatment of the subject can be found in Gantmacher 's two volume texts [12J from which 
we will here only cite the main definitions and results. For a more modern treatment, see Ref. 
|14| . Given two complex m x n matrices R and S, we form the homogeneous matrix polynomial 
^P(R,S) = M-R + i n variables \i and A. Two pencils "P(_r,s) an d ^(r^s') are strictly equivalent if 
there exists invertible matrices B and C independent of fi and A such that fiR' + AS" = B(nR + 
AS)C T . It immediately follows that V(r,s) an d 'P(r' j s') are strictly equivalent if and only if there 
exists invertible B and C such that BRC T = R' and BSC T = S'. 

The rank of Vmg\ is the largest r such that there exists an r-minor of V(r^) n °t identically zero 
(not equaling zero upon any complex substitution for \i and A). For % < r we let -D«(/x, A) denote the 
greatest common divisor of every i -minor of Vms) which is monic with respect to A. The invariant 
polynomials of pencil V(r,s) are the homogeneous polynomials Ei(/j,, A) = jjr^^l ^ or * = ^••■ r 
where Do(/j,, A) = 1. There will be a unique factorization of D r (^i, A) as D r (n, A) = /j, r ~ k pi ■ ... • pk 
where pj is of the form /j,Xj + A for xi G C, and the invariant polynomials will likewise have a 
factorization in terms of the pi and powers of /i. Then for each distinct pi (suppose there are q < k 
of them) we can consider the ordered set {p^ 1 p\ ir } where tij gives the largest power of pi that 
divides invariant polynomial Ej. The multiset generated by letting i range from 1 to q is called 
the elementary divisors of V^g)- Likewise, the set {/U e 'i, fi e ' r } is called the infinite elementary 
divisors of V^s) where e[ is the largest power of \i that divides Ei(n, A). From knowing the 
infinite and elementary divisors of Vr^s)/ its invariant polynomials can be constructed and vice 
versa. 

The right null space of Vr,s is the set of homogeneous polynomial vectors 

Xi(/x, A) = } j x i jff i ~ j X j 

j=0 
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such that (fj,R + A5)Xj(/z, A) = and / 0. A basis for the right null space whose elements 
have degrees e\ < ... < e p is called fundamental if any other basis whose elements have degrees 
e'l < ••• < £p implies e« < e' { for all i. A important property of any Xi(fi, A) = Ylf=o x ij^ tl ~^ ^ 
belonging to a fundamental set is that the x% are linearly independent fll2|. Likewise, we can define 
the left null space to be homogeneous polynomial vectors satisfying V] l s Xi(fi, A) = and form 
fundamental sets having degrees v\ < ... < v q . The values e\,...,e p and v\, u q are called the 
minimal indices of Vr,s- hi particular, the number of 6j that are zero will be called the zero index 
number, and the number of v\ that are zero will be called the transpose zero index number. With 
this overview, we can now state the main theorem characterizing strictly equivalent pencils. 

Lemma 1 (Kronecker). Two matrix pencils are strictly equivalent if and only if they have the same ele- 
mentary divisors (finite and infinite) and the same minimal indices. Moreover, suppose V(R t s) has finite ele- 
mentary divisors {(fj,xi + A) 611 , (fxxx + A) 612 , (/ux g + X) eqr }, infinite elementary divisors {n e>1 , ...,/i e ' r }, 
minimal indices ei,...,e p and u\, u q , a zero index number of g, and a transpose zero index number ofh. 
Then Vms) zs strictly equivalent to the canonical block-form diagonal pencil 

{0 h ^,L eg+1 ,...,L ep ,Ll h+i ,...,Ll,J} (1) 

e+l 

, * V 

/A fj, ... 0\ 

where hxg is the h x g zero matrix, L e = ^ ^ ^ "' ^ |e , and J full rank square pencil with 

\0 ... A iij 

block-form 

{N e '\...,N e '-,M eil ,...,M e " r } (2) 

where N e 'i = //F« + \H e 'i and M e ^ = (fixi + A)F« + fj,H e ^ with I* the t x t identity matrix and H l a 
t x t matrix whose only nonzero elements are ones on the super diagonal. 

We close this section by noting that both the minimal indices and elementary divisors of a pencil 
can be determined by efficient algorithms. The first involves determining the null space of scalar 
matrices, and the latter amounts to performing Gaussian elimination on the matrix fiR + XS [12\. 



3 Connection to 2 (g> m ® n Pure States 

Any 2® m®n state can be expressed inbra-ket form as | Vl/) = \0)a\R)bc + \^)a\S)bc- By choosing 
local bases {\i) s}i=o...m-i and {\i)c}i=o...n-i f° r Bob and Charlie respectively, we can express the 
state as 

|tf > = (\0) A (R ® I) + \1) A (S ® TjJ |*„) = (\0}a(I ® + ® S T )) |* m ) (3) 

where i?^ = a^, <Sjj = /%, and |<J>/A = Yli=o K)sK)c- Thus, there is a one-to-one correspondence 
between a2gm0n pure state |\&) and the pair of matrices (R, S), so that to every |^) and choice 
of indeterminates /i, A, we can uniquely associate the pencil Vms) which we shall equivalently 
denote as Vy. 
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There exists a nice relationship between the structure of and the local ranks of each sub- 
system. The reduced states of Bob and Charlie are obtained by performing a partial trace on the 
matrix |^)(^|. From above, then, it follows that 

PB = tr AC mm) =RR ] + SS^ 

Pc = tr AB m(*\) = R T R + S T S. (4) 

Here, "t" denotes the matrix transpose with respect to the basis \i)bc{]\ and "-" the complex 
conjugate of its entries. Also note that since Alice has a two dimensional system, her subsystem 
will either have full rank or be completely separated from Bob and Charlie. Combining these facts, 
we can prove the following. 

Lemma 2. (i) Bob and Charlie share pure entanglement (Alice separated) if and only ifVq can be expressed 
as a matrix polynomial in one indeterminate A; i.e. 

pR + XS = XS, (5) 

and (ii) Bob and Charlie's local ranks arem — h and n — g respectively where g is the zero index number of 
and h its transpose zero index number. 

Proof, (i) Alice is unentangled if and only if up to an overall phase, the state can be written as 

|0)(|i?) + \S)) + a\l)(\R) + \S)) which happens if and only if its associated pencil is (p + \a)R + 
(p + Xa)S = (p + Xa)(R + S) = XS. (ii) By definition, the zero index number is the number of 
linearly independent constant vectors \vi) such that R\v{) = S\vi) = 0. In this case, we must also 
have R\vi) = S\vi) = 0. It follows from Q that pc\vi) = if and only if R\vi) = S\vi) = and since 
complex conjugation does not affect linear dependence, we have rank(pc) = n — g. An analogous 
argument shows that rank(ps) = m — h. □ 

We now want to observe the effect of local invertible operators implemented by Alice, Bob, 
and Charlie; i.e. an SLOCC transformation. Any such operation can be decomposed as (A 
Ibc)(^-A ®B®C) where Bob and Charlie first act, and then Alice follows alone. When Bob and 
Charlie perform the invertible operator B (g> C, it is easy to check that the transformation | R) B c 
B <S) C\R) bc corresponds to R — > BRC T and likewise for S. Thus, the action of Bob and Charlie 
initiates the matrix pencil transformation pR+XS — > B(pR+XS)C T . In other words, local invertible 
operators of Bob and Charlie map matrix pencils to strictly equivalent ones. 

Any invertible operation by Alice can be represented by a matrix ( ° b d ) with ad — bc^ 0. Then 
the most general action by Alice will transform the state |^) as: 

\0)a\R)bc + \1)a\S)bc - \0} A (a\R) B c + c\S) BC ) + \1)a(J>\R)bc + d\S} BC ). (6) 

Hence, the corresponding pencil transformation is pR + XS — > (pa + Xb)R+ (pc + Xd)S = pR + XS 
where p = pa + Xb and A = pc + Xd. 

What concerns us is how the transformation (p, A) — > (p, A) affects the elementary divisors 
and minimal indices of a given pencil. For the latter, care must be taken since minimal indices 
are defined by the degree of polynomials in variables p and A. Nevertheless, the following lemma 
shows minimal indices to be an SLOCC invariant in 2 ® m <g> n systems. 
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Lemma 3. The minimal indices of a given pencil remain invariant under the action of Alice. 

Proof. Under an invertible transformation (/i, A) — * (jl, A) = (a/j, + bX, cp, + dX), a polynomial r- 
component vector p(p, A) = Ya^q 1 Sj=o x iil^ 1 ^ is identically zero iff p(p, A) = 0. To see this, we 
can introduce the standard basis {ek}k=i...r and consider p(p, A) as an rmre-component vector in 
the space spanned by basis p l X^e^. Then the transformation (p,, A) — > (jl, A) induces a homomor- 
phism on this space which thus cannot map any nonzero zero vector to zero. Consequently for any 
set of polynomial vectors {x.;(/i, X)}i=i... n (a) (p,R + XS)x.i(p, A) = iff (p,R + XS)xi(fa, A) = 0, and 
(b) {xj(/i, A)}i=i... n is linearly independent iff {xj(/t, A)}j=i... n is linearly independent, where linear 
independence means that for polynomials {p%(p, A)}j=i... n , Y^h=iP^^ ^) x «(M) A) = =>■ Pi(p, X) = 
for all i. Next, we claim that (c) for any set of linearly independent scalar vectors {x{j}j=o_ e with 
Xi € / 0, the highest degree of A having a nonzero vector coefficient in Xj(/i, A) = X^=o Xijtf ^ 
is the same as that in x;(/i, A) = J2'j=o Xi j^ t ~^' > '• ^ n ^ s f° nows because the coefficient of A e in 
Hj=o Xijfi € ~^XP = X]j=o x ij( a l J - + ^A) e ~ J (c/i + dXy is Yl e j=o x ijb € ^^d^ which is non-vanishing due 
to the linear independence of {xij}j=o... e - 

From (a), (c) and the linear independence of {xij}j=o...ei noted in the introductory discussion 
for any fundamental set of vectors, xi(ju, A) is a minimum degree polynomial in the null space of 
/j,R + XS iff xi(/l, A) is a minimum degree polynomial in the null space of fiR + XS. Now suppose 
that {xj(/i, A)}j=i... n are the first n vectors in a fundamental set for ijlR + XS iff {xj(/i, A)}j=i... n are 
the first n vectors in a fundamental set for fiR + XS. Then by (c), fj,R + XS and jlR + XS will have 
the same first n minimal indices. From (a), (b) and (c) again, x n+ i(/i, A) will be the next vector in 
the same fundamental set for fiR + A5 iff x n+ i(/x, A) is likewise for fiR + AS". Hence by induction 
and by running the exact same argument on (/iR + AS*) T , the lemma is proven. 

□ 

As for the elementary divisors, the situation is more complex since Alice's transformation can 
induce a mixing between infinite and finite divisors. By direct substitution, it follows immediately 
that after normalization, the divisors transform as 

, r ( „; + Ar; if mo 

\fi e i if 6 = 0, \fj, e » iib Xi + d = 0. 

We see that depending on the choice of A, infinite divisors can become finite and finite can become 
infinite. More importantly, given any general state having finite elementary divisors {(^uxj + A) 6 ^}, 
it is always possible for Alice to perform an invertible operation such that a / and {bxi + d ^ 0} 
for all i. As a result, we see that 

Any matrix pencil is SLOCC equivalent to one having no infinite divisors. (8) 

This observation simplifies the following analysis considerably since the general problem of deter- 
mining general SLOCC equivalence is reduced to the problem of equivalence among states having 
only finite elementary divisors. Combining the previous observations with Lemma [lj we arrive at 
the following theorem and a main result of this article. 
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Theorem 1. Two 2 <g) m <g> n states \ip) and \<p) having only finite elementary divisors {(fixi + A) e ^} and 
{(yi + A)A?} respectively are SLOCC equivalent if and only if their corresponding pencils are of the same 
rank, have the same minimal indices, = fijfor all i, j, and there exists a linear fractional transformation 
(LFT) relating the x% and m; i.e. for all i 



axi + c 
bxi + d 



Hi (ad-bc^O). 



(9) 
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yi 


1 






X\ 


Vl 




X\ 
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X2V2 
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Vl 
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X2V2 


x 2 


1 




X3V3 


2/3 
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X3V3 


X3 






X3 


2/3 


1 




X3V3 


x 3 


1 



A nice property of LFTsis that given any two trios {xi,X2,x 3 } and {2/1,2/2, 1/3} each with distinct 
values, there always exists a unique LFT relating the sets HJ. The form of the transformation is 
given by the determinants 



(10) 



We are now able to present an algorithm for determining whether two general 2 ® m <g> n pure 
states |?/>) and \4>) are SLOCC equivalent. 

(I) Input pencils and and determine their rank, minimal indices and elementary divi- 
sors. As noted above, this step can be achieved via polynomial-time algorithms. If the rank 
or minimal indices are not the same, \t/j) and \<f>) are inequivalent. Otherwise, perform an 
arbitrary LFT on them so the pencils only have finite elementary divisors {(/iXj + A) e,J } and 
{p,yi + respectively. By observation (|8) this can always be done. 

(II) Fix any three distinct x-i corresponding to divisors of powers e^. Choose any sequence of 
three distinct yi whose corresponding powers satisfy = eij and determine the LFT relat- 



ing (xi,yi), {x2,y2), and (£3,2/3) according to ( |10) . Choose a new Xi and determine if the 
LFT relates it to any remaining m belonging to an elementary divisor of the same power. By 
uniqueness of the LFT, if there is no such yi, the states are not equivalent. If there is, choose 
another x\ and repeat the search on the remaining y^ 

(III) If a perfect matching exists for all Xi and yi, then the states are equivalent. If not, repeat step 
(II) by choosing another ordered trio of the yi. If no LFT exists for all possible trios, the states 
are not equivalent. 

The Kronecker canonical form of an m x n pencil can be computed in time 0(m 2 n) (see the 
algorithm by Beelen and Van Dooren 111). F° r se ts of t elementary divisors, Step (II) this algorithm 
will require at most 0(i 3 ) steps. Thus the total running time is 0(m 2 n+mm{m, re} 3 ). Furthermore, 
the algorithm is constructive in nature because if two states are SLOCC equivalent, we determine 
the the specific a,b,c,d constituting Alice's operator in the transformation \cp) <slocc \ib). The 
operators Bob and Charlie are to perform can be determined from the invertible matrices that bring 
pencils and to their canonical forms of ([l} and are so-obtained by a Gaussian elimination 
procedure ||T2|. Hence, not only does our algorithm determine whether two states are equivalent, 
but it provides the necessary operators achieving the transformations. 



8 



4 All Tripartite Systems with a Finite SLOCC Equivalence Partitioning 



To count and characterize all the orbits, we will essentially find what combination of minimal in- 
dices and elementary divisors fit in an m x n matrix of form Q. A few simplifications will assist 
in this process. First, since any m x n pencil is simply the matrix transpose of an n x m one, it is 
enough to just consider m <n. Next, for a given dimension, we must only study the equivalence 
classes with Bob and Charlie having maximal local ranks since any rank deficient case will corre- 
spond to a class of maximum local ranks in a smaller dimension. To this end, Theorem fallows us 
to immediately determine the local ranks associated with each equivalence class. Furthermore, as 
evident from the Schmidt decomposition of any state with respect to bipartition AB:C, Charlie's 
local rank cannot exceed the product of Alice and Bob's. Consequently, if n > 2m, any state of a 
2 (8) m ® n system is the same as one in a 2 ® m ® 2m system up to a local change of basis on Charlie's 
part. This means that for the task of finite enumeration, we only need to consider systems up to 
dimensions 2 (g) 2 <g) 4 and 2 <8> 3 (g> 6. 

One further property of each equivalence class that we are able to study is the tensor rank. 
The tensor rank of a state is the minimum number of product states whose linear span contains 
the state, and this quantity turns out to be invariant under invertible SLOCC transformations 
|TlO| . For bipartite systems, the tensor rank is equivalent to the Schmidt rank, and a non-increase 
in Schmidt rank is also a sufficient condition for SLOCC convertibility between two such states; 
SLOCC equivalence classes are characterized completely by the Schmidt rank. Interestingly, in 
three qubit systems, tensor rank is also sufficient to distinguish between the various equivalence 
classes. However, we find that even for systems having a finite partitioning, the tensor rank is 
an insufficient measure for determining SLOCC equivalence. Our results follow from previous 
research on the tensor rank of matrix pencils done by Ja' Ja' fTBf and rederived in Ref. |5|. 

Lemma 4. [18], Let V(r,s) be a pencil with no infinite divisors in canonical form ^ with minimal 
indices e±, e p and v\, v q and J an I x I- sized pencil. Furthermore, let 5(J) denote the number of 
invariant polynomials containing at least one nonlinear elementary divisor. Then the tensor rank ofV^s) 
is given by 

X> + i) + X>j + 1 ) + z + ^ j )- ( n ) 

i=l j=l 

A summary of all the equivalence classes is provided in Table[T]in Appendix. We see that there 
are 26 distinct SLOCC classes for 2 ® 3 ® n (n > 6) systems. This reproduces the findings of Chen 
et al. [6] here obtained in an entirely different way by using matrix pencil analysis. 

5 Non-Invertible Transformations 

A natural question is whether it is possible to transform from one class to another via non- 
invertible transformations. One obvious constraint is that states with full local ranks cannot pre- 
serve their ranks under a non-invertible transformation. Consequently, we cannot convert among 
the states belonging to the same dimensional grouping above. A possible conjecture might be 
that unidirectional convertibility is achievable if none of the local ranks increase and at least one 
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decreases; certainly three qubit systems satisfy this hypothesis. This, however, is false in general 
as we will now observe. 

Let be some state having maximal local ranks of (2, m, n) and suppose \(j>) is a state with 
ranks (2,m,n — 1). If \4>) <slocc Alice and Bob's matrices inducing the transformation will 
be full rank while Charlie's will have rank n — 1. As for the latter, any such operator can be 
decomposed into a series of elementary column operations on followed by a mapping of the 
n th co l umn t a linear combination of the first n — 1; this n — 1-columned subpencil corresponds 
to the target state | <fi) . ^From the commutation relations of elementary operations, if we neglect 
permutations, 1 0} will be some n — 1-column subset of the original pencil following the application 
of just column-multiplying and column-addition transformations. Moreover, if column i is the 
linearly dependent column, then immediately after all column-additions of the i th column are 
performed, the remaining n — 1 columns must be equivalent to \<p). As a result, we obtain the 
following criterion. 

Theorem 2. Let and \(f>) be states with local ranks (2,m,n) and (2,m, n — 1), and let ci,...,c n 
denote the columns of V^(fi, A). Then \<p) < S locc \i>) ff for some 1 < i < n, there exists constants 
at, .., Oj-i, di+i, a n an d some invertible linear transformation (//, A) — > (/t, A) such that the pencil 
^Vi(A) A) = [ci + aia, c n + a n Ci] is equivalent to V<p([i, A). 

In general, for transformations in which Charlie's rank decreases to n—k, one need only modify 
this theorem by considering subpencils of V$ having n — k columns where to each of the columns 
is added a linear combination of the k non-included columns. Likewise, to account for transfor- 
mations when Bob's local rank decreases, the above criterion can be applied with the analysis 
conducted on the rows of V,/,(fJ>, A) instead of its columns. 

On the surface, Thm. |2]has limited value since it involves a search for values ai,..,ai-i, a,i+i,...,a n 
over the complex numbers. However, in many cases, it is easy to see whether or not such a collec- 
tion of numbers can be found. For example, for 1 < i < 4 in (ABC-19), upon any choice of the a,j 
and any transformation p,, A, the resultant pencil V^^fi, A) will either be rank two or it will have 
an elementary divisor of degree at least one. However, the state (ABC-18) is rank three with no 
non-trivial elementary divisors. Thus, the transformation (ABC-19)^(ABC-18) is impossible. On 
the other hand, for the state (ABC-17), when i = 1, we have detP^O, A) = A [A 2 - fi(^X + ^/x)] 
for 03 / 0. The state (ABC-8) has detV<p(ti, A) = A(/i+A)(2/i+ A). By choosing C2 = | and C3 = —\, 
these polynomials become equal as well as the elementary divisors, the ranks, and the minimal 
indices of the pencils. Thus (ABC-17)^(ABC-8) is achievable by SLOCC. 

In a manner similar to that just described, we have used Thm. |2]to analyze all possible trans- 
formations among the 2 ® 3 % n equivalence classes. Figure ([T]| in Appendix depicts the SLOCC 
hierarchy among the classes. 

6 Conclusions and Future Research 

In this article, we have used the theory of matrix pencils to study 2 <8> m n pure quantum states. 
In doing so, we were able to derive a polynomial time algorithm for deciding SLOCC equivalence 
of such states. For all tripartite systems having a finite number of equivalence classes, we have 
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obtained state representatives and determined the partial ordering among these classes based on a 
criterion for general SLOCC convertibility in 2 ® m (g) n systems. It is interesting to note that in the 
hierarchy chart of Fig. 2, there exists certain transformations that are impossible even though the 
local rank of Charlie decreases by two. The transformation (ABC-14) to (ABC-7) is such example. 

A natural extension of this work is to find efficient algorithms for deciding LOCC equivalence, 
LOCC convertibility and SLOCC convertibility in 2 (g m (g n systems. We have made progress on 
those questions. Another natural next line of inquiry might to consider p(g m<g n systems and their 
corresponding degree p matrix polynomials. Indeed, much research has been conducted on higher 
degree elements, especially those having special properties such as being symmetric |T4|. Unfor- 
tunately, there exists no corresponding characterization like Kronecker's for strict equivalence of 
matrix pencils of degree greater than two. Making the project of generalizing to higher degrees 
more dubious is the fact that determining SLOCC equivalence for pgmgn can be reduced from 
a tensor rank calculation on a set of p bilinear forms [7j, and this problem has no known solution 
for p > 2 [5 1 (the general problem is, in fact, NP-Hard ||T6| ). 

As noted in the introduction, we are not the first to study SLOCC convertibility in multipar- 
tite systems, and it would be interesting to try and develop the relationship between our results 
and the work of others. For example, Miyake's results involve "hyperdeterminants" and their 
singularities [21 j. It would be valuable to investigate the correspondence between matrix pencils 
and hyperdeterminants or to introduce the connection to the quantum information community 
if such a correspondence has already been obtained. In another work, Liang et al. have recently 
proven a set of conditions both necessary and sufficient for the convertibility of two qubit mixed 
bell-diagonal states [19J. As these mixed states can be considered pure with respect to a 2 (g 2 ig 4 
system, it would be fruitful to study transformations between tripartite "purified" bell diago- 
nal states via our matrix pencil construction and compare it to the convertibility conditions in 
Ref. [19|. Doing so might suggest ways in which purified tripartite pencils can assist in deciding 
equivalence between general 2 <g n mixed states. 
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Appendix A 



Representative 


Local Ranks 


Tensor Rank 


Representative 


Local Ranks 


Tensor Rank 


(A:B:C) 


(144) 


1 


(AB:C) 


(2,2,1) 


2 


(AC:B) 


(2,1,2) 


2 


(A:BC-1) 


(1,2,2) 


2 


(ABC-1) 


(2,2,2) 


2 


(ABC-2) 


(2,2,2) 


3 


(ABC-3) 


(2,2,3) 


3 


(ABC-4) 


(2,2,3) 


3 


(ABC-5) 


(2,2,4) 


4 


(ABC-6) 


(2,3,2) 


3 


(ABC-7) 


(2,3,2) 


3 


(ABC-8) 


(2,3,3) 


3 


(ABC-9) 


(2,3,3) 


3 


(A:BC-2) 


(1,3,3) 


2 


(ABC-10) 


(2,3,3) 


4 


(ABC-1 1) 


(2,3,3) 


4 


(ABC-12) 


(2,3,3) 


4 


(ABC-13) 


(2,3,3) 


4 


(ABC-14) 


(2,3,4) 


4 


(ABC-15) 


(2,3,4) 


4 


(ABC-16) 


(2,3,4) 


5 


(ABC-17) 


(2,3,4) 


4 


(ABC-18) 


(2,3,4) 


4 


(ABC-19) 


(2,3,5) 


5 


(ABC-20) 


(2,3,5) 


5 


(ABC-21) 


(2,3,5) 


6 



Table 1: Equivalence Classes in 2 ® 3 <8> 6 Systems 
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2®3®6 




2®2®4 



2®3®3 



aim: 17, 

ABC 11, 

ABC-8, 
ASC-2 



ABC 12, ABC 11, 
ABC-10, ABC-9, 
ABC 3, ABC 2 



ABC 13, ABC 12, 
AIM! 11, ABC-10, 
ASC-2 



ABC-13, aim: 11, 
ABC 9, ABC 3, 
ASC-2 



ABC-10, 
ABC-9, 
ABC 2 



ABC-12, 
ABC 11, 
ABC 8 



ABC-13, ABC-12, 
ABC-11, ABC-10, 
ABC-9, ABC-8 



ABC-12, 
ABC 11, 
ABC 8 



ABC-13, ABC-12, 
ABC-11, ABC-10, 
ABC-9, ABC-8 




2lg>3®2 



Figure 1: Complete Hierarchy of SLOCC Equivalence Classes; Arrows Indicate a Non-invertible 
Transformation. 
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